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SUMMARY 

This paper illustrates how the use of a general analysis gjJJg. 

and analymg fluid/mechanical systems. One such package » 

Swvices product The basic transmission line equations for moUmg 1 ? S^vasY5^ 

SS ^paper describes how to analysis tool has been used to model several 
fluid subsystems of the Space Shuttle Orbiter. 

INTRODUCTION 

Modeling complex fluid/mechanical systems can involve difficulties ^yond d^bmg 
the syarn «3y. N« only docs An mk involve coding of An Ktud nqMMns, An 
analyst is also facnd wiA numnrical integration of Aosn cquaoonsjUsOTO^on trf An 
Lid rat-processing of An results. Thus, Ante ousts a need for a toolwhrch 
SSnL. into a singln package 

analvsis program has been found to be one such tool which can be used to ^ectiveiy 
mffiSSStanical systems. With the advent of fast workstations based on RISC 
Sto. graphically interfaced analysis programs for system analysis are highly efficient. 

Modeling using EASY5 can be done in a finite-element type 
subroutines. The user defines the behavior of a single element * 1 ^- 

pipe flow or a spring-mass system) using the appropriate 
discretizes the system as a combination of these elements, 

metod routines The features of to code benefit the user by no^wr and 

linear analysis capability. Nonlinear time-domain simulations cante^ 
several different integration methods. This package also has the "^ ***"$* 
system to provide transfer function, root locus, eigenvalue, as wdfl f 

analysis. Also contained within EASY5 is a plotting routine which can provide plots of 

results for the different types of analysis. 

While any of the systems that could be modeled using EASY5 
using FORTRAN, this type of software represents a convenient combination of many of 
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the tools which the analyst requires and significantly reduces time required to develop a 
new system simulation. 
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Figure 1. Flow chan of EASY5 modeling process. 

THEORY 

Most fluidAnechanical systems can be classified into subsets of similar components: 
pipe flow, pipe intersections (tees and crosses), orifices, volumes, and spring-mass 
systems. If the modeler has subroutines defining these components, they can be combined 
to represent complicated systems. 

The basic building blocks for the fluid flow subroutines, or macros in EASY5 
terminology, are the transmission line equations (ref. 1). The equations are listed below 
(see Figure 2 for notation). 


Ijmj = Pj - P i+1 - Rfilmilmi (1) 

Cj Pj = mj_i - mj ; i = 1, N (2) 

where: 

If = inertance of the fluid element, 

Cj = capacitance of the i^ 1 fluid element, 

mj = mass flow into the i+1 element, 

Pj = pressure at the center of the i^ 1 fluid element 
Rfj = resistance, 

N = total number of fluid elements used to model a line segment 

For a uniform line modeled with equal-length elements, die inertance, capacitance, 
flow resistance and temperature equations are tire same for all elements and are given by 
(assuming one-dimensional flow and isentrqpic behavior): 
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where: 

L = fluid element length 
A = flow area 

Y = polytropic process exponent 
T = temperature 
R = gas constant 
f = friction factor (pipe flow) 

D = line internal diameter 
t = equivalent length for minor losses 

Equations 3, 4, and 3 
can also be expressed for j 
fluid. 



Ci 

Figure 2. Typical discretization of a fluid line segment 


specify the flow parameters for gas systems. These parameters 
i liquid system by using the bulk modulus and density of the 


The standard transmission line equations can be modified to 
and crosses by using additional flow equations. Volumes of changing size can be mode 
using Equation 7, which assumes an isentropic process. 


P« 7 (mRT-PVyV (7) 

The mass flow rate through an orifice is given by the familiar relationship (ref. 2) 
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where Cp = orifice discharge coefficient, A = orifice area, and the subscripts 1 and 2 

denote pres sure/temperature upstream and downstream, respectively. A valve can be 
modeled using a variable area orifice. To approximate the opening of the valve, it has been 
found th at varying the area using a hyperbolic tangent function yields die best results. 
However, any type of continuous or discrete function could be used as long as die rates of 
change within the model do not become too large for the integration step size. 

Because EAS YS requires systems of first order differential equations, spring-mass 
systems are modeled by breaking the system's second order differential equations into first 
order equations. For example, die governing differential equation far a spring-mass- 
damper system, 

x = -(cx + kx - F(t))/m (9) 

may be replaced by the following two first-order equations: 

v = -(cv + kx - F(t)Vm (10) 

X = V (11) 


APPLICATION 

An EASY5 macro is very similar to a FORTRAN subroutine. The macro contains the 
code required to describe the behavior of a single model element, e.g., a transmission line 
element, spring-mass combination, etc. The parameters which define the physical 
characteristics of the element are inputs to the macro, as are the boundary conditions for 
that element as calculated by an adjacent element. The outputs of the macro are the values 
calculated using the code within the macro and the specified inputs. A model is then built 
by linking a series of macros together using their inputs and outputs. 

For example, consider the three element section of a model shown in Figure 3. An 
acoustic line is being modeled using a macro named *TR' (EASY5 macro names consist of 
2 char acters) The acoustic line macro is a combination of the pressure/flow differential 
equations, isentropic temperature relationship, and a curve fit of the Moody diagram. The 
macro first calculates the current temperature assuming an isentropic process. Next, the 
macro uses a logic block to determine which way flow is moving. After die flow direction 
is determined, the friction factor is calculated using die Reynold's Number and the 
equations describing the Moody diagram. The flow and pressure derivatives are then 
calculated and integrated. These outputs are then used as inputs to other elements. 
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Figure 3. Acoustic Line Connection and a Simple Three Element Model 


When modeling an acoustic line, the line is broken up into a senes TJj c 

teneAof each element is based on the highest frequency of interest ^n^tf *© a 

lin^The usual FEM rales regarding the minimum and maximum number of elements m a 

pressure of the next line element must be inputs to the current line element 

A long length of line can require an excessive number of transmission 
In order to minSdze the effort required to build the model, a 
Uncmaao was developed. Thecode internal to 

user specifies the number of sub-elements to be model 

1 to 999? This development greatly reduces the amount of time required to develop 

of a system. 

The time step used for nonlinear time-dcmrin timutoiwt 
nature of the model. Hie optimum time step it found thi ou|h an ■^^go ^ tforfljted 
time step integrators, while variable time step integration schemes hareto^ra^jMnng 
totterarion tunTstep to maintain solution accuracy with the largest acwpable nmestep. 
x5 STcrfTtofsip results in numeric* emirdue toUr*e w 

ofa rimVstep can cause excess round-off enm. He opnmiun ti me f 

solutions has been found to be one which, when reduced, gives results 
SfSTswS. The recurrence formula foe the wave equations 
wtenchoosing a time step size. Therefore, the following relationship needs to be 

considered. 
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where At is the time step size. Ax is the element length and c is the speed of sound of the 
media being modeled. A detailed explanation can be found in reference 3. 

EAS Y5 offers several different types of integration algorithms. These include: 
Fixed-Step YmaMcSttP 

Euler (1st order) BCSGear 

Huen (2nd order) Adams-Moulton 

Fixed-Step Runge-Kutta (4th order) Stiff Gear 

Variable-Step Runge-Kutta 

The variable step integration schemes adjust the integration step size based on how fast 
the system states are changing. Ideally, these methods would be desirable for use since 
they represent a potential execution time savings. However, it has been the authors' 
experience that the variable step methods are not particularly compatible with the macros 
that have been developed to model fluid systems, due to die quadratic damping term and the 
large pressure derivatives associated with small elements. Typically, the integrator ends up 
iterating excessively trying to optimize the step-size, thereby greatly increasing the 
execution time. Good results have been obtained using the variable-step methods on 
spring-mass systems. 

Another nonlinear analysis feature of EASY5 is steady state analysis. The steady state 
command returns the equilibrium operating condition of die model, the model rates of 
change are essentially zero for this analysis. 

EAS Y5 is also capable of linear dynamics analysis. This is done by linearizing the 
state equations in the model by pertuibating them about the operating point to create a linear 
perturbation model. This linear model can then be used for other types of analysis such as 
transfer function, root locus, closed loop eigenvalue and other frequency domain analyses. 


EXAMPLE 1 - 750 PSIA MPS HELIUM SUPPLY REGULATOR 




In this example the authors were asked to investigate a problem with the Space Shuttle 
main propulsion system (MPS) 750-psia helium pressure regulator. Two regulators 
experienced full-open failures due to high frequency (900 Hz), b’gh amplitude oscillations. 
The failures took place on a new test stand which was constructed to replace die original 
regulator qualification stand after it was destroyed in the collapse of the building in which it 
was located. 

The authors were tasked to develop dynamic models of die test stand as well as models 
of all three MPS engine helium supply system configurations utilizing an existing model of 
the regulator developed by the vendor. The purpose of the models was to determine the 
source of the oscillations, evaluate potential for oscillations on the Orbiter, and to test 
possible solutions for correcting the problem. 
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Moling Effort 

UK EASY5 software was seated for this be 

La the regulator were modeled by using bi-linear springs. 
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Figure 4. MPS helium supply regulator schematic 


Analysis 


Using the transfer function option of EASY5, it was determined thu 

ISSTan? pressure oscillations occurring downstream of the regular m this fiequency 
range. 

A model of the complete newly constructed verification test stand 
oscillatory behavior of the regulator was duplicated using the 
mSched the first acoustic mode downstream of 

confieuration's fundamental frequency coincidentally matched that of the regulators 
bdtow^^^h leadto fatiguefailure of the bellows. Models of the complete Orbiter MPS 

helium supply system were also constructed (Figure 5). Each engme supp y y 

consisted of approximately 1000 degecs-of “freedom. 
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PANEL B PANEL A 



Figure 5. Diagram of Engine 1 helium supply system 


The EASYS model of the regulator was used to guide and evaluate design changes 
proposed by the vendor. The final design showed stable operation in both tests of actual 
hardware and in numerical time simulations with the math model. Figures 6 and 7 show a 
Bode plot and simulation results of the regulator before and after the redesign. 

■ Original Design • New Design 




Frequency (Hz) 

Figure 6. Bode plot of regulator outlet pressure to controller pressure transfer function. 
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Figure 7. Time simulations of the regulator before and after redesign. 


EXAMPLE 2 - PRCS THRUSTER 

This second example describes the use of the EASY5 software in 
propellant feedline dynamic models of the Orbits ScmodeS Solved 

g^s^ waLrtammer mulysis «» . complex two-ptase flow of the chug 

stability of the thruster. 

Background 

Combustion stability testing of the PRCS thruster involves injecting helium into the 
propellant feedlines in order to provide a combustiM disturban^ 
anrerisely known quantity. The test stand is designed to produce a nominal Aw of . 
hKdlmng steady stare conditions. However due to 

the flow of helium into the thruster can va^ vndely over nme^ For this “ 

model was desired to predict the amount of hehum mgesred bythe jhroster. f which 

propellants for the thruster are monomethylhydrazme and nitrogen tetroxide, both of whic 
are liquids at the operating pressure and temperature. 


Modeling 

The test stand models were developed using macros similar to those used for the MPS 
heUa^Sp^jS The models ^comprised of single- md two-ptase etanems. 
line elements severi inches upstream of die helium injection point arc capable of two- 
phase flow representation, while the remainder of the transmission line elements are singl 
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phase (see Figure 8). The two-phase macros assume a homogeneous gas-liquid mixture, 
ideal and isentropic behavior of the gas phase, and are based on equation (13): 



YmjP 
■yRTm- + 

8 Pi 




(13) 


where 

B = liquid bulk modulus 
Vg ■ gas volume of element 

mg = mass flow of gas 

mj = mass flow of liquid 
Pl = density of liquid 



Figure 8. Schematic of PRCS Test Stand Model 

Initial simulations used a time history table of thruster dumber pressure measured in 
test firings as the boundary condition at the end of the propellant feedline. Based on the 
steady state flow rate and pressure drop, the resistance of the line could be fine-tuned to 
achieve the required flow parameters. Time-domain simulations used fourth-order Runge- 
Kutta as the integration method, with an integration step size of 1.0E-0S seconds. 
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Results from these models were in the form of time histOTy ^loo 
parameters. Sample results are shown in Figures 

behavior of fuel feedline pressure just upstrcamof*e thrasttr. J^^nnjBxaoofiran^hMe 
Snniarions was also used to size the test stand fcedline lengths 10 obtain a wiiertiainmer 
fteauencv similar to the actual vehicle installation. Figure 10 shows thebebum mgestjon 

160 millisecond 

such thatthe ending conditions for one simulation are Ae £jtial “ndi^toften^ 

Sr*. PRCS modeling task, this allowed a sequence of firings to be simulat ed so thatth e 

Sount of heUum resirfing in the lines would build as the test firing sequence progressed. 



Figure 9. Typical s i m ulatio n result. 
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Figure 10. Helium injection profile. 


fnr thei^ps regulator task Variable area orifices were used to represent the opening and 

Sg va^l^a^as the valve poppet moved. The stiffness of *e poppet spring is 
reDreS ented by tabular data taken from tests conducted during the valve tevetopment 
program. Leak rates around the poppet seals are simulated by not allowing the variable 
area orifices to close completely. 

A diagram of the model schematic is shown in Figure 11. Test stand ^branon, which 
may contribute to some of the high amplitude pressure and acceleration oscillations 
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observed during tests, is included in the model. The test stand is treated as a single degree- 
of-frcedom system having mass, stiffness, and damping characteristics close to that of the 
stand. Test stand motion is applied to the fluid system through macros which have moving 
boundaries. The test stand velocity is applied to these elements as an element wall velocity, 
which drives the element pressure derivative. Also included in the model is the combustion 
timelag. The flow out of the last injector element is delayed from combusting (expanding 
into gas) by a specified amount of time. This is accomplished through an EASY5 
continuous delay macro. The chamber pressure is calculated based on the capa ci ta n ce of 
the chamber, the amount of fuel and oxidizer flowing into the chamber, and the amount of 
gas flowing out of the chamber. The amount of gas flowing out of the chamber is 
determined using the characteristic velocity (c*). 

Due to the small size of the injector, very small elements were necessary to obtain the 
required fidelity. The size of these elements dictated that the integration step size also be 
small. The optimum step size was found to be IE-07 seconds, using fourth-order Runge- 
Kutta as die integration method. 



n» chanter prawn obtained Irom tha oombtabon moM warm m 
boundary oendtton lor thainjaalore. Tha tvurt • UMd aa • forcing function lor 
tha M aland maaafelfntaa modal. 


Figure 1 1 . Schematic of detailed PRCS thruster model. 
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CONCLUSIONS 


^ss^sszssssssss^^ 

nro^or has proven to be effective and convenient for modeling toese types ( of wtems. 
E£f^ litv models of several complex non-linear fluid, structural, and mechanical y*ten 
S^dSns were developed which correlated weU wi* test data and provided a basis for 
analyzing and eliminating causes of adverse dynamic interactions. 

NASA- JSC Propulsion Branch is continuing to use EASY5 for other propulsion 
svstemsA subsSd S“f macros and models have been developed wbch«Uwquick 
^^uiate Sytical results to be obtained for a wide variety of propulsion fluid and 

mechanical systems. 
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